PHOTOVOLTAIC AIRCONDITIONING, 
REFRIGERATION AND HEATING; AN 
ENVIRONMENTAL FRIENDLY OPTION 


by 

Kirti Kumar Dhawan 





DEPARTMENT OF ELECTRICAL ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY KANPUR 

September, 1999 


PHOTOVOLTAIC AIRCONDITIONING, 
refrigeration and HEATING: AN 
ENVIRONMENTAL FRIENDLY 

OPTION 


A Thesis Submitted 

in Partial Fulfillment of the Requirements 
for the Degree of 
Doctor of Philosophy 


by 

Kirti Kumar Dhawan 


to the 


DEPARTMENT OF ELECTRICAL ENGINEERING 
INDIAN INSTITUTE OF TECHNOLOGY KANPUR 

September, 1999 






Certificate 


Certificate 

It is to certify tliat the work contained in the thesis entitled Photo- Voltaic Aircon- 
ditioning, Refrigeration and Heating: An Environmental Friendly Option 
by Mr. Kirti Kumar Diiawan, has been carried out under my supervision and that this 
work has not t)een submitted elsewheie for a degree. 


Dr. Manohar Prasad 
Professor 

Department of Mechanical Engineering 
IIT Kanpur 


S(ipteml)er, 1999 


1 



Acknowledgements 


It is a rare opportunity to have Professor Manohar Prasad as my thesis supervisor. I 
am highly indebted to him for the inspiring guidance and useful discussions I had with 
him throughout the course of this work. 

I am highly indebted to Professor G.K.Dubey for his guidance, suggestions and 
technical feedback in doing the literature survey and finally suggesting this problem 
and his advice for solving the same. 

I am also grateful to Dr. B.G. Fernandes for his guidance and help in purchasing 
photovoltaic array and batteries. His time to time comments and suggestions has 
immensely luilped in solving this problem. 

My sincere thanks to Dr. U.S. Dixit, Asst. Prof., Mechanical Engineering Depart- 
ment, I.I.T. Guwahati for his sincere and consistent help throughout this work. 

Thanks are due to Shri P. P. Bhola for his timely help and full cooperation in setting 
up and conducting the experiment. 

I wish to thank Shri R.K.Bajpai for his excellent drawings/ tracings and Shri A.K.Sriv: 
for his error free typing. 


Kirti Kumar Dhawan 



Publications 


Following publications have been communicated out of this work: 

Accepted 

Following publications have been accepted: 

1. ‘Application of Fuzzy Set Theory to Photovoltaic Airconditioning’ in 1999 In- 
ternational Appliance Technical Conference, Florida, U.S.A. Under Consideration for 
publication in IEEE Transactions: Industry Applications. 

2. ‘Photovoltaic Airconditioning: An Environmental Friendly Option’ in 20th In- 
ternational Congress of Refrigeration, IIR/IIF, Sydney, 1999. 

Under Review 

Following publications are under review and consideration for publication: 

1. ‘Photovoltaic Heating’ in Journal of Solar Energy Engineering, Milwaukee, 
Wisconsin, U.S.A. 

2. ‘Photovoltaic Airconditioning and Heating: Optimum Load Matching of a 
Battery Connected Photovoltaic Array’ in 2000 International Appliance Technical 
Conference, Florida, U.S.A. 



3. ‘Photovoltaic Refrigeration’ in Indian Journal of Engineering and Materials 
Sciences, National Institute of Science Communication, New Delhi, India. 


IV 



Synopsis 


The sources of conventional energy are dwindling fast with a corresponding rise in cost. 
Therefore, considerable attention is being paid to other alternative sources of energy. 
Solar energy is one of the non-conventional energy resources which has the advantage 
of being freely available in plenty, particularly in India. It can be harnessed as thermal 
energy for water heating, cooking, absorption-refrigeration and solar thermal power 
generation. It can be converted directly into electricity by using photovoltaic array. 
A photovoltaic array, called PV array, consists of a number of solar cells arranged in 
series and parallel. The electrical energy generated in the array can be stored in a 
battery and can be used to run a compressor of a vapor-compression system or to heat 
the generator of a vapour absorption system. 

The airconditioning, refrigeration and heating require energy in the form of heat 
(in the generator of a vapour absorption system) or in the form of work (to run the 
compressor of a vapor compression system). They require mostly the same type of 
appliances and process and differ in application. Whereas, in the context of hot coun- 
tries like India, airconditioning is understood as a process of cooling a surrounding, 
the heating means the reverse of it in which a space is heated by the rejected heat of 
condenser. The refrigeration is defined as the artificial withdrawal of heat, producing 
in a substance or within a space a temperature lower than that which would exist under 
the natural influence of surroundings. Thus, refrigerating machine is an integral part 
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of airconditioning system. 

The term photovoltaic airconditioning, heating and refrigeration means running 
these systems by the solar energy converted in the form of electricity. As the solar 
energy is not available round the clock and also the weather conditions are highly 
uncertain, the electrical energy is stored in the battery. Apart from the continuous 
depletion of natural resources and degradation of environment due to burning of fossil 
fuels, the solar energy is very much relevant for remote area applications. However, at 
present, the solar energy is very costly and therefore, the size of the array should be 
selected judicially. 

In the present work, a computer code is developed to find out the array size de- 
pending on the application. A model has been proposed for finding out the array 
photo-generated current. The model is validated by conducting experiments on a 100 
W photovoltaic array at IIT Kanpur. The cooling load, heating load and refrigeration 
load are computed on timely basis based on ASHRE handbook. Feasibility of two types 
of systems- vapour compression and vapour absorption system is considered. It is seen 
that the vapour compression system is economical. In this, the electrical energy will 
be used to run a D. C. motor which will be used to drive the compressor. 

As the weather conditions are highly uncertain, the fuzzy set theory is employed 
to find out the array size. In this, the weather parameters are expressed as linear 
triangular fuzzy numbers and possibility distribution for array size is obtained. Then, a 
method has been proposed to find out the array size based on the subjective requirement 
of the designer. 

Once a proper size of the array is chosen, it has to be matched with the load 
connected to it. The load is a battery of 12 V or 24 V. By load matching is meant 
the proper selection of series and parallel rows of cells for a given total number of 
cells. A load matching factor is defined as the ratio of energy stored over the day to 
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the maximum possible energy stored over a day. The number of series and parallel 
rows of cells and internal resistance of the battery is found by solving the optimization 
problem to maximize the stored energy with the constraint of total number of cells and 
generated voltage being greater than the battery voltage. A parametric study has been 
carried out and it is seen that a load matching factor of more than 0.95 may easily be 
obtained. However, the system will be highly inefficient in case of a load mismatch. 

The study has been also carried out to find out the proper inclination of the array 
in order to get maximum output from the insolation. For this purpose also, an opti- 
mization algorithm has been employed. It has been also suggested to cover the space 
by photovoltaic array in a manner in which it will fulfill the twin objective of shading 
the roof to reduce the cooling load (refrigeration requirement) and generating electricity. 
Finally, the environmental benefits of using photovoltaic systems are studied in detail. 
It is seen that use of solar energy will reduce carbon dioxide level and ash particulates 
in the environment. 

Thus, the present work studies various aspects of the photovoltaic airconditioning, 
refrigeration and heating and discusses the importance of its use. The use of a new 
model for array-photo-generated current, application of the fuzzy set theory and a 
methodology to find out the optimum load matching and inclination of the array and 
environmental study are the main features of the present work. 
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Chapter 1 

INTRODUCTION AND 
LITERATURE REVIEW 

1.1 Introduction 

Airconditioning is defined as the simultaneous control of temperature, humidity, clean- 
liness and air motion. Depending upon the requirement, air conditioning is divided into 
the summer air conditioning and the winter air conditioning. In some hot countries of 
the world, the term air conditioning is usually employed for summer air conditioning 
and for winter air conditioning the term heating is used. 

The refrigeration may be defined as the artificial withdrawal of heat, producing in 
a substance or within a space a temperature lower than that which would exist under 
the natural influence of surroundings. A refrigerating machine is an integral part of 
the air conditioning system, because if the surrounding temperature is higher than the 
room temperature, heat is to be removed from the space to be air conditioned. How- 
ever, for the places where surroundings are at a temperature lower than the required 
temperature, the heating has to be done. The machinery which performs the heating 


1 



process is called a heat pump. The main difference between the refrigeration system 
and heat pump can be physically conceived of from the fact that in the former there is 
pumping of heat out of the system as against pumping of heat from surroundings into 
the system in the latter. Thus, a refrigeration system can be used as a heat pump by 
suitable arrangements, which would cause refrigerant to flow in such a way that heat 
is transferred to the conditioned space and is drawn from the surroundings. 

The vapour compression refrigeration system is the most widely used refrigeration 
system. It basically requires electrical energy to run a compressor. The compressor 
raises the pressure of the refrigerant vapour, thus raising its boiling point. The vapour 
is condensed in a condenser. The condensed liquid is throttled through a valve or a 
capillary tube into the evaporator. As a result, the pressure drops causing the boiling 
point of the refrigerant to drop. Therefore, the refrigerant evaporates (and takes heat) 
at much lower temperature thus cooling and dehumidification of the conditioned space. 
The low-pressure vapour is sucked by the compressor and cycle continues. 

Figure (1.1) shows the schematic diagram of a refrigeration system along with 
process of energy supply to the system. In order to produce electricity to run the 
compressor, the fuel is burnt in the boiler. The steam coming from the boiler enters a 
steam turbine which produces mechanical energy to run the generator. The generator 
produces electricity which is transmitted to various places. There are losses at every 
stage. 

In refrigeration, an important term called refrigeration effect is defined as the 
amount of cooling produced by a system. This cooling is obtained at an expense 
of some form of energy. Hence, it is customary to define a term known as Coefficient 
of Performance (COP): 

Refrigeration effect 
COP — ' r 7 

energy input 
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Fig. 1.1 Schematic diagram of a refrigeration 
system « 







For a heat pump the term. PI (Performance Index), is defined as, 

pj _ Heating effect ^ 
energ}- input 

If the same system is used as a refrigerator and a heat pump with operating conditions 
unchanged, then PI of the heat pump will be COP of refrigerator plus one. 

In the case of vapour-compression system driven by electricity, COP is based on 
electrical energy input. However, when comparing two systems employing different 
forms of energies, total energy input (which includes losses) should be taken into ac- 
count. For example, in the arrangement of Figure (1-1), COP calculated on the basis 
of energy supply through fuel may be just about one third of that calculated on the 
basis of energy supply through electric input to compressor. 

The vapour-absorption refrigeration system eliminates the need for a compressor. 
Figure (1.2) shows a basic absorption unit. Low-pressure vapour from the evaporator 
is absorbed by the liquid solution in the absorber. If the absorption process were 
executed adiabatically, the temperature of the solution would rise and eventually the 
absorption of vapour would cease. To perpetuate the absorption process the absorber 
is cooled by water or air that ultimately rejects this heat to the atmosphere. The pump 
receives low-pressure liquid from the absorber, elevates the pressure of the liquid, and 
delivers the liquid to the generator. In the generator, heat from a high-temperature 
source drives off the vapor that had been absorbed by the solution. The liquid solution 
returns to the absorber through a throttling valve whose purpose is to maintain the 
pressure difference between the generator and absorber. The Coefficient of Performance 
of the absorption cycle COP is defined as. 


+ COP 


COP = 


Refrigeration rate 
Rate of heat addition to generator 


(1.3) 


The heat input to the refrigeration system can be given by various means including 


solar energy. 
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Highpressure vepor 



FIG. 1.2 THE BASIC ABSORPTION UNITS 




Fig, 1.3 100 W Photovoltaic Refrigerator 






It is clear from, the above description that the air-conditioning, refrigeration and 
heating systems essentially require energy in the form of heat or work. As the sources 
of conventional energy are dwindling fast with a corresponding rise in cost, considerable 
attention is being paid to other alternative sources of energy. The solar energy, which 
is one of the non-conventional energy resources, has the advantage of being pollution 
free and being freely available in plenty, particularly in India. It can be converted into 
thermal energy or directly into electricity (by photovoltaic conversion) . A photovoltaic 
array, called PV array, consists of a number of solar cells arranged in series and parallel. 
The electrical energy generated in the array can be stored in a battery and can be used 
to run a compressor of vapour compression system or to heat the generator of vapour- 
absorption system. Figure (1.3) shows the block diagram of the refrigeration system 
using PV array. The PV array is a very costly item, therefore its size should be selected 
judicially in order to meet the load requirement in the presence of uncertain weather 
conditions. Also, the effort should be to utilize the PV array most effectively. 

The present work investigates the various aspects of photovoltaic air-conditioning, 
refrigeration and heating. Photovoltaic systems are not being used on a commercial 
scale due to extremely high cost. It is expected that in the coming century they are 
expected to attract more attention in view of the global warming and degradation of 
environment as a result of burning of fossil fuel and depletion of fossil fuel once and 
for all. The present work aims at providing guidelines for the design of systems using 
PV array. 


1.2 Literature Review 

There is hardly any literature available on the use of PV array in air-conditioning, 
refrigeration and heating systems. However, matching of PV array with the loads and 



photovoltaic, in general, has been studied by various authors. The present section gives 
a brief account of the research carried out in this area. 

Optimum matching of directly coupled loads to PV array is studied [Khouza- 
m,1990]. Matching of resistive loads coupled directly to PV array is optimized by 
selecting the load rated power, voltage and resistance with respect to those of the ar- 
ray [Khouzam,1990]. But direct matching of other loads such as D.C. motor has not 
been studied. Optimum matching for separately excited D.C. motors can be obtained 
by changing of connections of cells and modules within the solar cell array [Braunstein 
and Zinger, 1981]. This can only give discrete matching. Two to three different series - 
parallel arrangements of modules can give up to 15% increase in efficiency. More than 
three different states (connections) are not economically viable [Zinger and Braun- 
stein, 1981] and [Salameh and Dagher,1990]. Continuous optimum matching can be 
obtained by changing module connections along with field current variation [Faldella 
et ah, 1991]. 

Choppers are also used for continuous matching purpose. A simple step-up chopper 
circuit can consist of a single power transistor and an inductor. There exists a unique 
duty ratio which gives optimum utilization efficiency at all insolation levels. This 
property makes this technique attractive since it eliminates the continuous adjustment 
of the duty ratio ■'092]. Optimum matching can also be obtained by 

predetermining the flux constant K of the motor for a particular load type and the 
location of the installation before placing the order for the same. Due to the losses in the 
motor, the maximum gross mechanical power characteristic is different from maximum 
power characteristics of the array. Flux constant K is determined [Saied,1988] to 
operate the motor at its maximum mechanical power points for different insolations 
and then K is weight-averaged to adjust for different duration of insolation intervals. 
D.C. motor with this K can be directly coupled to PV array for optimum matching. 
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Similar study is carried on in [Saied,1991] for a D.C. motor connected to PV array 
via a chopper. Here output energy is plotted with respect to flux constant K and the 
optimum range for K and corresponding value for duty ratio of the chopper is decided. 
Comparison of operation between various types of D.C. motors shows that separately 
excited D.C. motors are most suitable for such systems. Efficiency of D.C. motors for 
part load or at light load is high. Also since the source is D.C., therefore D.C. motors 
become the obvious choice to be run by PV array. But brushes in D.C. motor need 
replacement (after 4000 h of operation) otherwise commutator may get damaged. 

A D.C. motor can be replaced by (Inverter + Induction motor) with improved part 
load efficiency arrangement. In reference [Bhat et ah, 1987] induction motor/pump 
system is connected to PV array via a chopper and an inverter (VSI, Voltage Source 
Inverter). Here maximum power is extracted for any insolation condition and optimum 
V/ f relationship is used for the motor control which is not quite correct. 

As mentioned earlier, maximum mechanical output power characteristic of the mo- 
tor is different from the maximum electrical power characteristic of the array and also 
optimum V/f relation will be different as pump head (load) varies. Induction motor- 
pump system can also be connected via a chopper and CSI (Current Source Inverter) 
[OJO,1991] instead of VSI. 

Reference [Enslim and Snyman,1991] suggests that traditional cascaded PV power 
conditioning system, consisting of many stages such as maximum power point tracker, 
regulator, battery and an inverter can be replaced by a composite circuit comprising all 
the above functions in a circuit. But the circuit given here is only for constant frequency 
output whereas for improving part load-efficiency of induction motor, we need variable 
frequency supply. Proper selection of solar array and battery size, known as sizing 
problem [Jaboori et ah, 1991] is based on annual simulation using actual meteorological 
(insolation and temperature) and typical load data pertinent to the place of installation. 
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The procedure determines the useful (effective), dumped (wasted) and commercial 
energy components in the system. Optimum PV- and battery sizes are selected where 
effective energy is maximum. 

The residential photovoltaic power conditioning system is studied in [Bose et al.,1985]. 
High frequency link scheme is used here, giving considerable weight reduction of the 
power converter and maintaining a unity power factor at the AC line terminal. Mi- 
cro computer is used for array current feedback control, maximum power tracking, 
sequencing and safe zone steering control. 

[Ranade et al.,1989], present islanding problem - a phenomenon in which one or 
more inverters continue to operate and keep the system energized even after the utility 
supply is disconnected. It is shown that islanding can be virtually eliminated through 
the use of destabilized feedback controls for the inverters. 


1.3 Scope and Objective of the Present Work 

As is evident from the review of literature, the load matching of the array has not 
been studied for air-conditioning, refrigeration and heating. The present work aims at 
finding the array size for specified design parameters. A code has been developed to 
calculate load of the system as a function of time. The size of the array is selected 
seeing the load requirement in a day. As the load as well as the available PV array 
output keep changing with time, it is advisable to always employ a battery. Thus, load 
matching problem becomes a problem of a battery connected with PV array. 

The present work also incorporates the uncertainties in the weather conditions using 
fuzzy set theory. A methodology for designing the array size using fuzzy set theory is 
suggested. This approach will be helpful to produce an economical and reliable design. 

For a given size of the array, the cells can be arranged in the series and parallel in 
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various ways. They should be arranged in a manner to give maximum output over the 
day. One objective of the work is to find out the optimum arrangement of the cells 
in the PV array. The PV array should be positioned in a manner such that it gives 
maximum output and covers the roof properly reducing the load due to radiation. This 
work addresses these aspects. 

Finally, the gain to environment due to use of PV array operated systems is dis- 
cussed. At present the solar energy may not be competitive with the conventional 
sources except at remote places where other sources of energy can reach only with 
difficulty. However, if the pollution aspect is considered, the solar energy may be a 
good alternative to conventional sources of energy. 

It is shown that if in place of conventional electricity photovoltaic energy is used 
then for a room size as given in Appendix G, CO 2 at the rate of 3.75 kg/s at the peak 
time can be saved. Also for the same size of room amount of ash particulates saved is 
0.8 kg/h. 


1.4 Organization of Thesis 

The thesis is organized in the following manner. Chapter 1 presents introduction, liter- 
ature review and the objective of the work. The methodology to find out the array size 
is presented in Chapter 2. The methodology to find out cooling/heating/refrigeration 
load as a function of time is discussed here. Chapter 3 is on the application of fuzzy 
set theory for estimation of array size. Chapter 4 is on the optimum arrangement of 
series and parallel cells in the array. The positioning of array and the environmental 
issues are discussed in Chapter 5. Chapter 6 is on conclusions and scope for further 
work. 
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Chapter 2 


DETERMINATION OF 
PHOTOVOLTAIC ARRAY SIZE 

2.1 Introduction 

Airconditioning, refrigeration and heating essentially require same appliances. In these, 
either a vapour-compression system or a vapour-absorption system is employed [Prasad, 
1993]. There is a requirement of electricity for running the compressor in the vapour- 
compression system and pump in the vapour-absorption system. The required electrical 
energy may be obtained by photovoltaic array. Since the photovoltaic array is a costly 
item, its size for a specific purpose has to be selected judicially. 

In this chapter methodology to find out the size of the array is described for aircon- 
ditioning, heating and refrigeration. Both, vapour-compression and vapour-absorption 
systems have been considered. A software has been developed for finding the total 
heat load requirement and the size of the array. The determination of the size of array 
is the most important part of the photovoltaic airconditioning, heating and refrigera- 
tion design as this gives the idea about the space occupied by the array and the cost 
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required. Subsequent chapters will be dealing with the arrangement of the array. 

2.2 Characteristics of Photovoltaic Array 

The photovoltaic array output varies with solar insolation and temperature. The array 
characteristics has been described in reference [Khouzam, 1990]. The V-I equation of 
a solar array generator is expressed by: 

V = AVr. ln{(Ip^ -1 + J,)//,} - LR, (2.1) 

where, I is the load current, Iph is the array photo-generated current, Ir is the total 
reverse saturation current, Rg is the total series resistance and AVp is the equivalent 
thermal voltage of the array. If the array generator is composed of Ng series cells and 
Np parallel strings, then. 


Ipk 

= h.Np 

(2.2) 

Ir 

= Is.Np 

(2.3) 

Rg 

= Rc.Ng/Np 

(2.4) 

AVt 

= Ng.Vp 

(2.5) 


where, Ii is the photo-generated current per cell, Ig is the reverse saturation current 
per cell, Rc is the series resistance per cell, and Vp is the cell thermal voltage. 

The effect of the array temperature on the reverse saturation current, the photocur- 
rent, and the thermal voltage is given by [Appelbaum, 1986] 


Ir{T) 

= 4(ro).(T/To)^exp{-6(l/T-l/ro)} 

(2.6) 

Ipk(T) 

= /p,(ro).[i + a(r-ro)] 

(2.7) 

AVt(T) 

= AVr(ro).r/To 

(2.8) 
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where, To is the cell reference temperature (at 288.15 K), and a and b are coefficients 
dependent on the type of semiconductor material used and the manufacturing process. 
The output power of the array depends on the loading and is given by 

P = LAVt. ln{IpH - / + Ir)/Ir} - I^.Rs (2.9) 

The array photo generated current varies with time over a day. Anis et al. [1985] 
and Appelbaum [1987] have taken a standard clear sky insolation model in which the 
following equation holds good; 

7’p/i = /g- sin(u)) (2.10) 


where, Iq is the array photo-generated current at solar noon, and w is the solar hour 
angle. Actually the variation of Iph over the day is not sinusoidal, rather it is propor- 
tional to the radiation falling on the array. In the present model, it is assumed that 
the following relation holds good for the Iph- 

(radiation at a particular time) 


Iph = Ig- 


( 2 , 11 ) 


(radiation at solar noon) 

The above model is validated by doing the experiments on an inclined photovoltaic 
array. Figures (2.1-2) compare the power output obtained experimentally and its com- 
parison with the theoretical prediction using the proposed model. It is seen that the 
power prediction from the model is in good agreement with the experimental results. 


The next section describes the procedure for finding the solar radiation at a partic- 
ular time. The V-I and P-V characteristics of PV array at different level of insolation 
are shown in Figure (2.3). As I is varied, V varies, therefore for a particular (Vmp, Imp) 
power input is maximum for a given insolation. The array maximum power is found by 
differentiating power (P) in equation (2.9) with respect to load current I and equating 
to zero. This gives the following equation: 





Time ( hour ) 

Fig. 2.2 PV array output, Experimental and theoretical (little cloudy and 
dusty day in May) 


Power (watt) ^ 

N>u>.psaic»^oou> 






Array current (A) 
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Fig. 2. 3 Effect of Insolation on the I-V and P-V 
characteristics of crystalline silicon 
solar cells. 


Array generated power(W) 





Iph — Imp + -^r{exp[(2/TOpi?5/ylVT) + Imp! (Iph — Imp + -^r)] “ 1} (2.12) 

The maximum power point voltage, Vmp can then be found from equation (2.1) by 
replacing I by Imp- The maximum power, Pmp is then given by: 


P —VI 

mp ^mp'-^mp 


(2.13) 


For the efficient utilization of PV array, the array should be operated at maximum 
power point voltage. 


2.3 Radiation Falling on the Surface 

Based on an analysis of U.S. data, ASHRAE [1982] has given a method for estimating 
the hourly variation of global and diffused solar radiation falling on a horizontal surface 
on a clear day. The relationship may be developed which find the radiation on an 
inclined surface. It is advantageous to describe some definitions first. If 6 is the angle 
between an incident beam of flux / 6 „ and the normal to a plane surface, then the 
equivalent flux falling normal to the surface is given by Ibn cos 6. The angle 6 can be 
related by a general equation to (j) the latitude, S the declination, 7 the surface azimuth 
angle, u the hour angle, and /? the slope. Each of these will first be defined. 

The latitude </> of a location is the angle made by the radial line joining the location 
to the centre of the earth with the projection of the line on the equatorial plane. By 
convention, the latitude is measured as positive for the northern hemisphere. 

The declination 5 is the angle made by the line joining the centers of the sun and 
the earth with its projection on the equatorial plane. It arises by virtue of the fact that 
the earth rotates about an axis which makes an angle of approximately 66.5° with the 
plane of its rotation around the sun. The declination angle varies from a maximum 
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value of + 23.45° on June 21 to a minimum value of -23.45° on December 21. It is zero 
on the two equinox days of March 21 and September 22. Cooper [1969] has given the 
following simple relation for calculating the declination: 


5 (in degrees) = 23.45 sin 


- 360 

.365 


(284 -h n) 


(2.14) 


where n is the day of the year. 

The surface azimuth angle 7 is the angle made in the horizontal plane between the 
line due south and the projection of the normal to the surface on the horizontal plane. 
By convention, the angle is taken to be positive if the normal is east of south and 
negative if west of south. 

The hour angle w is an angular measure of time and is equivalent to 15° per hour. 
It is measured from noon based on local apparent time (LAT), being positive in the 
morning and negative in the afternoon. The term ‘local apparent time’ is defined below. 


Local Apparent Time - This time can be obtained from the standard time ob- 
served on a clock by applying two corrections. The first correction arises because of the 
difference between the longitude of a location and the meridian on which the standard 
time is based, the correction has a magnitude of 4 minutes for every degree of difference 
in longitude. The second correction called the equation of time correction is due to the 
fact that the earth’s orbit and rate of rotation are subject to small perturbation. This 
correction is based on experimental observation and is plotted in Figure (2.4). Thus, 

Local apparent time = Standard time =F 
4 (Standard time longitude — longitude of location) 

-[-(Equation of time correction) (2-15) 

The negative sign in the first correction is applicable for the eastern hemisphere while 
the positive sign is applicable for the western hemisphere. 
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Fig. 2.4 Equation of time correction 



Another way of defining solar hour angle is like this, 


w = 15(t — 6) 


( 2 . 16 ) 


where w is in degrees, t is the solar time, which is related to actual time ta by: 


t = 6.0 + 


6.0{ta - tsr) 
(j'ss tgr') 


(2.17) 


where and tgs are the sunrise and sunset times respectively. The sunrise and sunset 
times at a particular place are given by: 


tsr = 12 — cos ^ (— tan (?!>. tan J) (2-18) 

15 

tss = 12 + -^cos~’^(— tan^i. tanJ) (2.19) 

15 

Here 0 is the latitude of the place and 5 is declination, given by equation (2.14). 

The slope /3 is the angle made by the plane surface with the horizontal. It is taken 
to be positive for surfaces sloping towards the south and negative for surfaces sloping 
towards the north. 

The relationship between various angles is : 


cos 6 = sin 0(sin 5. cos /? + cos 5. cos 7 cos w sin p) 

+ cos 4){cos 6 cos u; cos /5 — sin 6. cos 7 sin P) 

+ cos 5 sin 7 sin w sin ,5 (2.20) 


For horizontal surface, P is zero and 

cos 6 = sin (/>. sin 6 + cos (p. cos 6. cos w (2.21) 

The angle 6 in this case is called the zenith angle and is denoted by Oz. 

The global radiation Ig reaching a horizontal surface on the earth is given by 

Ig=Ib + Id . (2.22) 
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where 


Ig = hourly global radiation 
h = hourly beam radiation 
Id — hourly diffused radiation 
Now, 

h — IbnCosB^ (2.23) 

where Itn is the beam radiation in the direction of the rays and is the angle of 
incidence on a horizontal surface, i.e., the zenith angle. Thus, 

Ig = Ibn cos 9^ + Id. (2.24) 

In the ASHRAE model, it is postulated that for a clear cloudless day 

Ibn = Aexp[— R/ cos^z] (2.25) 


and 

Id = CIbn (2.26) 

where A, B and C are constants whose values have been determined from an analysis of 
data for the U.S.A. and are given in [Sukhatme, 1994], and are reproduced in Appendix 

A. 

Radiation on a Tilted Surface: The radiation on a tilted surface is composed 
of three parts, beam radiation, diffused radiation and the reflected radiation. 


(1) Beam radiation: The ratio of the beam radiation flux falling on a tilted surface to 
that falling on a horizontal surface is called the tilt factor Rb for beam radiation 
and is given by: 


Rb = 


cos 9 
cos9z 


(2.27) 
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(2) Diffused Radiation: The tilt factor for diffused radiation is the ratio of the 
diffused radiation flux falling on the tilted surface to that falling on a horizontal 
surface. It is given by: 


Rd = 


1 + cos ^ 
2 


(2.28) 


It is to be noted that (l+cos/?)/2 is the radiation shape factor for a tilted surface 
with respect to the sky. 


(3) Reflected Radiation: Assuming that the reflection of the beam and diffused radi- 
ations falling on the ground is diffused and isotropic, and that the reflectivity is 
p, the tilt factor, Rr, for reflected radiation is given by 


Rr = 


P 


T — cos^' 
. 2 . 


(2.29) 


The flux It falling on a tilted surface at any instant is thus given by 


It — IbRb + IdRd + {Ib + Id)Hr (2.30) 

2.4 Estimation of Cooling/Heating/Refrigeration Load 

The design of an airconditioning/heating/refrigeration plant requires corresponding 
load to be estimated with sufficient accuracy. The load calculation needs the informa- 
tion about climate and inside conditions desired. A generalised computer programme 
for hourly load has been developed. The procedure is similar to that given in ASHRE 
[1982]. 

For the ambient temperature variation, a mathematical model [Alford et ah, 1939] 
has been adopted which gives a close approximation to an average value equal to the 
mean daily temperature range. Therefore, the final form of the expression comes out 
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to be a harmonic series of first order. 

rp f'^min Tma,x\ fT-max Tmin\ rir/'j. j. t\i /n 01 \ 

To = j + j cos[l5(t - tsr - I)] (2.31) 

where To is the outside temperature and Tmin and T^ax are the minimum and max- 
imum temperature in a day. In deriving the above relation, the assumption is that 
minimum temperature over the day occurs one hour before the sunrise, and maximum 
temperature occurs 12 hours after that. This equation has been verified for Kanpur 
conditions by K. Sekar [1986]. 

The Sol-air temperature For heat transfer calculations, it is convenient to combine 
the effects of outside air temperature and solar radiation clustered into a single fictitious 
quantity called sol-air temperature (Tj). The rate of heat transfer, q, from outside to 
inside surface of a sunlit structure, considering the effect of emissivity (e) of the surface 
may be written as, 

Q — hoiTo — Tyj^o) + cifh ~ cAi? (2.32) 

where, 

ho is outside air convective heat transfer coefficient (W/m^ K), 

To is outside air temperature (K), 

Tyj^o is temperature of the outside surface(K), 

at is absorptivity of the structure, and It is total intensity of solar radiation (IT/m^K) 
€, emissivity of the surface, AR = 63 W/rri^ for the horizontal surface and zero for the 
vertical surface. The heat transfer can also be written as 

q = ho{Ts — Txo^o) (2.33) 

Comparing the equation (2.32) and equation (2.33), the expression for sol-air temper- 
ature, considering emissivity of the surface [ASHRE, 1982] is: 


Tg = To Ot-It/ho — cAR/Hq 


(2.34) 



When emissivity of the surface is considered, is found to be about 0.3 K less than 
that of the case when emissivity is ignored. The outside air convection heat transfer 
coefficient, hg is calculated as follows [Kadambi and Hatchinsen,1968]: 


hjQ 


14.278 H- 3.559V;; 
18.423 + 3.81K; 
26.55 + 5.066314 
28.64 + 6.364K; 


for very smooth surface 
for smooth surface 
for rough surface 
for very rough surface 


(2.35) 


where 14 is the design velocity of the outside air in km/h and hg is kW/m^ K. 

Selection of Indoor Conditions: - The field studies [Malhotra, 1967] reveal that in 
India, for hot and humid climates, the comfort zone is between 295 K and 297.5 K 
effective temperature. 

Here the comfort condition is taken as follows: 

Tdb = 297K, Rff = 0.65, V — 0.8 m/min. 

The components of cooling load which are considered here are as follows: 

(a) Heat Transfer Through Walls and Roofs - It is calculated by time lag and decre- 
ment factor method [Arora,1993]. Combined effect of outside air temperature 
and incident solar radiation intensity is accounted, as sol-air temperature is used 
for all calculations. Thus, the structural heat load Qst including solar radiation 
is found to be 

5 5 

Qst = T.Qk = E 

fc=l K=l 

where Tgk is the mean outside sol-air temperature on the surface of the 
structure, \k is the corresponding decrement factor, r is the time lag and is 
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the product of overall heat transfer coefficient and surface area. Tsk{ta — t) is the 
sol-air temperature at (ta — r) time. is the inside design temperature of the 
room and ta is the actual time at which cooling load is being calculated. 

To obtain the overall heat transfer coefficient, U, is found as follows. It has been 
assumed that the walls and roofs have three layers of materials with thickness 
Xi and thermal conductivity ki of the layer of the structural material. In the 
present case, the values are taken from ASHRE[1982]. U is given by 


U = 


J. |_ ^ 

ho 



Here, hi is the inside convective heat transfer coefficient. 


(2.37) 


(b) Heat Transmission Through Glass Windows - Solar radiation, direct and diffused, 
incident upon a glass surface is, in parts, transmitted, reflected and absorbed. 
Thus r, r and a represent the respective fractions known as transmissivity, re- 
flectivity and absorptivity, then 

r -f r -f a = 1 (2.38) 


The heat gain of a space through glass then comprises; 

1. All the transmitted radiations 

2. A part of the absorbed radiation that travels to the room, and 

3. The heat transmitted due to the temperature difference between the outside 
and inside. 


The direct radiation enters the space only if the glass is receiving the direct rays 
of the sun. The diffused radiation enters the space even when the glass is not 
facing the sun. 
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It can be shown [Arora, 1993] that the heat gain of space is 


0 = (A.„ rolo + AtM) + A.,naoIo + AaM 

1 + f 

+UA{To - Tin) (2.39) 

where Agun is the glass area directly exposed to the sun, A is the total glass area. 
Subscript D and d denote the terms for direct and diffused radiations respectively. 
U, the overall heat transfer coefficient is given by 

1-11 

u~ To 

The values of transmissivity and absorptivity of ordinary glass are taken from 
[Arora, 1993] and are reproduced in Appendix B. 

(c) Infiltration Load - Infiltration is the leakage of outdoor air into a building through 

cracks and opening caused by pressure difference across the boundary surfaces. 
The exchange of air may lead to both heat and moisture gain for the space. The 
volume of infiltration is related in terms of the room volume and tabular values 
are available in [Prasad, 1993] and are reproduced in Appendix C. The expression 
for the infiltration load is: 

<3infil = NACHiha - ha)/24 kW (2.40) 

■^air 

where Vroom is the volume of the room, is the specific volume of air at 
outside temperature, Nach is the number of air changes/day, ha and hia are the 
enthalpies of outside and inside air, respectively. 

(d) Ventilation Load - Ventilation air requirement can be taken from tables given in 
Table A.7 of reference [Prasad, 1993], reproduced in Appendix D. The procedure 
for calculating ventilation load is the same as that of infiltration load. 


27 



(e) Internal Heat Gains - It comprises occupancy load, lighting load and appliances 
load. Occupancy load is calculated from values given in Table 19.1 of Arora 
[1993] and reproduced in Appendix E. Lighting and appliance load is calculated 
by adding total wattage of lights .and equipments. 

In case of refrigerators, internal heat gain consists of loads due to storing of fruits 
and vegetables and for making ice. The information needed for calculating this 
is as follows: 

(i) The average amount of commodities present in the refrigerator 

(ii) The average amount of commodities which is replaced at each hour. 

(iii) The per hour ice production. 

The heat of respiration of products is given in Table 19.3 of [Arora, 1993] and 
reproduced in Appendix F. 

The total load at a particular time is found by adding all the components of load 
at that time. Seeing the hourly load, a suitable size of the system may be selected, and 
components can be designed accordingly. Integrating load with respect to time using 
Simpson’s formula, the energy requirement for 24 hours period is estimated, which is 
needed for the estimation of PV array size. 

2.5 Estimation of Size of PV Array 

Current Imp corresponding to maximum power can be found from equation (2.12). 
The Newton-Raphson method fails to give reasonable results, because first derivative 
of the function with respect to current is greater than 1. Interval reducing method 
[Pike, 1986] is used to solve equation (2.12), where absolute value of the function is 
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minimized for I^p values between .75Ip/i and .95Ip/i. In this method the convergence is 
guaranteed. The maximum power, P^p can then be found from equation (2.13). The 
Pmp is found at different times in a day and is integrated with respect to time using 
Simpson’s formula to find out energy available in 24 hours. 

Using DO loops, energy available in a day is computed for each combination of Ns 
and Np. The combinations which give total electrical energy in a day close to the elec- 
trical load requirement for refrigeration (refrigeration load/coefficient of performance) 
are printed. A suitable combination of Ns and Np can then be selected to meet re- 
frigeration requirement. Knowing Ns and Np one can find the space needed for PV 
array. 


2.6 Results and Discussion 

The array size was estimated for airconditioning, heating and refrigeration. The fol- 
lowing subsections discuss typical results for the estimation of photovoltaic array. 

2.6.1 Results for air conditioning 

The programme was run for the PV array estimation for airconditioning a room in 
Kanpur. The particulars of room are given in Appendix G. Figure(2.5) shows variation 
of cooling load with time for 24 hours period for the month of May (month of maximum 
heat), that is why May has been taken for design conditions for cooling load calculation 
[Kadambi and Hutchinson, [1968]) . It is clearly seen that cooling load varies drastically 
with time. Figure (2.6) shows the electrical energy requirement for vapour-compression 
and vapour-absorption system. For calculating energy requirement, values of COP 
for compression and absorption system are taken as 3.84 and 1.2 respectively, from 
section 7.1.3 and Chapter 7 of reference [Prasad, 1993]. The graph of electrical energy 
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requirement versus time for vapour-compression system has been obtained by dividing 
cooling load by COP, while for the vapour-absorption system it is assumed that starting 
from 3 hours after sunrise upto 2 hours before sunset all the power required can be 
obtained by solar collectors and therefore electrical power requirement is zero. Rest of 
the time, electrical power required is obtained by dividing cooling load by COP. The 
total electrical energy requirement in a day is found from the area under the electrical 
power requirement curve. From Figure (2.6), we find that total energy required is 
15.91 kWh for the vapour-compression system and 23.80 kWh for vapour-absorption 
system. A 1730 W array in the case of vapour-compression system and 2630 W array 
in the case of vapour-absorption system is needed in order to meet electrical energy 
requirement. The powers available from 1730 W and 2630 W arrays are shown for the 
month of May in Figure (2.7). The PV cell data are given in Appendix H. It is found 
from the computer program that 320 series cells and 22 parallel strings will be required 
for 1730 W array. For 2630 W array 320 series cells and 33 parallel strings are needed. 
The size(area) of the array required is estimated referring [Tata BP Solar India Ltd. 
Product Catalogue]. The array size for 1730 W array is 12m x 12m. For 2630 W array, 
it is 12m X 20m. Figure (2.8) shows the cooling load at different times in the month 
of September and Figure (2.9) is electrical power needed in that month for the two 
systems. The powers available from the 1730 W and 2630 W arrays are shown for the 
month of September in Figure (2.10). It is seen that powers available in this month 
are more than the requirement. The excess power is stored in the battery. 

Table (2.1) shows total energy required and energy available from PV array from 
April to October, for both vapour-compression as well as vapour-absorption systems. 
The energy required for vapour-compression systems is much lower than that of vapour 
absorption system because COP of vapour-compression system is higher than that of 
vapour-absorption system. It is seen that when array is selected on the basis of the 


30 






Ele 


Vapour N 

absorption' \ 

' \ 

I \ 

I \ 


2.00 


1.00 


Vapour 

compression. 


O.OOl 


8 12 
Time (h) 


Fig. 2. 6 Electric power requirement 




Vapour absorption 








Elcctric{KW) 


35 



Fig, 2. 9 Electric power requirement for September. 








Table 2.1: Total Energies (in kWh) per day required and available in different months 


for airconditioning. 


Month 

Vapour Compression System 
Energy required PV energy 
available 

Vapour Absorption System 
Energy required PV energy 
available 

April 

10.44 

14.99 

15.19 

22.48 

May 

15.91 

15.90 

23.80 

23.81 

June 

12.71 

16.12 

18.89 

24.18 

July 

13.59 

15.83 

22.21 

23.74 

August 

11.47 

15.04 

20.95 

22.57 

September 

12.59 

14.00 

20.14 

21.00 

October 

4.93 

12.81 

4.75 

19.21 


month of May, in other months power available is more than the requirement. This 
takes care of the uncertain weather conditions in those months. 


2.6.2 Results for Heating 

Programme was run for the PV array estimation for heating a room in Kanpur. The 
particulars of room are given in Appendix G. The variation of heating load with time 
of 24 hour period for the month of January (the coldest month) is shown in Figure 
(2.11). It is clearly seen that heating load varies drastically with time. Figure (2.12) 
shows the electrical energy requirement for the vapour-compression system. For cal- 
culating energy requirement the values of COP for heating system is taken as 4.84. 
In section 7.1.3 and chapter 7 of reference [Prasad, 1993] a typical COP value of a 
practical airconditioning system based on vapour- compression system is given as 3.84. 
Performance Index(PI) of heat pump is one plus COP of airconditioner. That’s why 
the value 4.84 is taken here. The total electrical energy requirement in a day is found 
from the area under the electrical power requirement curve. It comes out to be 11.9 
kWh. A 1796 W array is needed in order to meet this requirement. The PV cell data 
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are given in Appendix H. The array size (area) for this array is about 160 m^. Power 
available from this array is shown in Figure (2.12). 

Figure (2.13) shows the electrical energy requirement for vapour-absorption system 
for heating. Based on the COP value of 1.2 for airconditioner of vapour absorption 
type, PI here has been taken as 2.2. It is also assumed that starting from 3 hours 
after sunrise upto 2 hours before sunset all the power required can be obtained by 
solar collectors and therefore electrical power requirement is zero. Rest of the time 
electrical power required is obtained by dividing cooling load by COP. the total energy 
requirement in a day is 23.9 kWh. A 3615 W PV array is found to meet this array 
requirement. The size(area) of the PV array is 325 m^ . The power output from such 
an array is shown in Figure (2.13). 

Figure (2.14) shows heat load for the month of November. It is observed that heat 
load is less compared to January. Consequently, the power requirements in vapour- 
compression and vapour-absorption systems are less (Figures (2.15-16)). The PV arrays 
have already been chosen based on the requirements in the month of January. Their 
outputs for the two systems are also shown in Figures(2. 15-16). It is seen that PV 
array output is not much different from the output for the month of January. That 
means when the array is designed on the basis of the coldest month, in other months 
the system is highly over designed. Table (2.2) clearly shows this. In this table energies 
required and ayailable for all months of winter are shown. It is seen that in photovoltaic 
heating, there is an inverse relation between required and available power which is 
obvious because increasing solar radiation increases PV array output but decreases 
heat load. Therefore it is not economical to design the system based on the coldest 
month requirement. 

In the places where all the year-round airconditioning is employed, the array size 
which is sufficient for cooling the space can be used for heating also. In the previous 
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Table 2.2: Total Energies (in kWh) per day required and available in different months 
for heating 

Month Vapour-Compression System Vapour-Absorption System 
Energy required PV energy Energy required PV energy 
available available 


January 

11.9 

11.9 

23.9 

24 

February 

11.3 

13.3 

22.3 

26.8 

March 

8.4 

14.0 

15.4 

28.1 

November 

6.7 

12.9 

14.5 

25.8 

December 

9.7 

12.6 

19.5 

25.4 


example, the array size was computed for the airconditioning of the same room in the 
hottest month. Interestingly, the array size for the vapour -compression system comes 
out to be very close to the array size selected here. Therefore, it is advisable to employ 
the same size of array for the airconditioning and heating purpose. In that case, in 
some months, array will be supplying more power than the requirement. Whenever 
the available power is more than the required power, the excess power is stored in a 
battery which takes care for the peak requirement periods. 


2.6.3 Results for Refrigeration 

The programme was run for the PV array estimation for a refrigerator of 45 liter 
capacity. The particulars of refrigerator are given in Appendix I. The variation of 
refrigeration load with time for 24 hour period for the month of May (month of max- 
imum heat) is shown in Figure (2.17). It is clearly seen that refrigeration load varies 
with time. Figure (2.18) shows the electrical energy requirement for the refrigerator. 
The refrigerant chosen here is R-12. For calculating energy requirement the values of 
COP for refrigeration system is taken as 3.5 as calculated in Appendix J. The graph of 
electrical energy requirement versus time for refrigerator has been obtained by dividing 
refrigeration load by COP. 
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The total electrical energy requirement in a day is found from the area under the 
electrical power requirement curve. From Figure (2.18), we find that total energy 
required is 897 Wh. A 100 W array is found to meet this energy requirement. The 
power output from this array is shown in Figure (2.19). This array gives 904 Wh in a 
day. The 40 series cells and 10 parallel strings will be required for this array and space 
required is 0.96m^ . The PV cell data are given in Appendix H. 

Figure (2.20) shows the refrigeration load at different times in the month of January 
(the coldest month) and Figure (2.21) is electrical power needed in that month. The 
power available from the 100 W is shown in Figure (2.22). It is seen that power available 
in this month is more than the requirement. The excess power is stored in the battery. 

Table (2.3) shows total energy required and energy available from PV array for all 
months. It is seen that when array is selected on the basis of the month of May, the 
same array is sufficient for meeting the energy requirement in other months also. The 
code was run for the environmental conditions of Bombay also. The electrical energy 
required in the month of May comes out to be 784 Wh there, while the array is able 
to supply 850 Wh in a day there. In the month of January energy required is 715 Wh 
and array supplies 716 Wh in a day. Thus the same 100 W array is able to supply 
the required power in Bombay also. The two cities selected here are the representative 
cities where environmental conditions and latitude differ much from each other. It is, 
therefore, expected that in other cities too the 100 W array is expected to be sufficient. 
The specifications of refrigerator is given in Appendix I. It is considered as maximum 
capacity specifications. The refrigerator is not used to its full capacity all the time. 
Moreover, the heat of respiration, frequency of door opening for replacing the food 
items etc. have been taken on higher side. In practice the energy requirement is 
expected to be less than calculated here. The excess energy will be stored in battery 
to be used during uncertain weather conditions. 
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Table 2.3: Energies required and available per day for 100 W refrigerator. 


Month 

Energy Required (Wh) 

PV Energy Available (Wh) 

January 

599 

663 

February 

622 

710 

March 

683 

771 

April 

835 

839 

May 

897 

904 

June 

864 

916 

July 

872 

900 

August 

852 

855 

September 

814 

816 

October 

725 

727 

November 

659 

672 

December 

621 

647 


2.7 Summary 

A method to estimate array size for airconditioning, heating and refrigeration system 
is presented and computer code has been developed for this purpose. For each of these 
cases, examples have been taken for estimating array sizes. For airconditioning and 
heating both vapour-absorption and vapour- compression systems were considered. It 
is observed that array size requirement is much more in the case of vapour-absorption 
system than in the case of vapour-compression system for both airconditioning and 
heating. This is due to the low COP of the vapour-absorption system. 

For the case of a refrigerator, 100 W, 45 liter capacity refrigerator can be developed 
which can successfully operate at Kanpur and Bombay, the two representative cities of 
India where environmental conditions differ much more from each other. 
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Chapter 3 


APPLICATION OF FUZZY SETS 
TO THE ARRAY SIZE 
ESTIMATION 

3.1 Introduction 

As discussed in Chapter 2, the array size is dependent on the cooling load required and 
the solar energy available. Environmental conditions constitute major uncertainties 
apart from other type of uncertainties, such as number of occupants inside the room. To 
account for uncertainties, we are familiar with two types of theories : probability theory 
and fuzzy set theory. Probability theory is used when an event is precisely defined and 
an adequate statistical data are available. Various parameters are tested as stochastic 
variables and it is assumed that uncertainty and imprecision are of random nature. 
However, most of the meteorological phenomena have another kind of uncertainty 
referred to fuzziness due to imprecision or lack of data and subjectivity of opinion 
or judgement. These are sources of uncertainty which are caused by vague definition 
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rather than by chance. Fuzzy sets deal with the membership or non-membership of an 
object in a set with imprecise boundaries. 

In meteorology, vague informations often arise, for example, the sentence ‘today is 
hot’ consists in fuzziness. In this case, there is no exact definition about which degree 
of temperature is considered as ‘hot’. For weather forecasting, especially for long range 
forecasting, we usually do not ask about the detail but instead the tendency of future 
weather. Therefore, it seems more reasonable to treat the array size determination 
problem from the fuzzy set theory [Zadeh, 1965; Kaufmann and Gupta, 1985; Klier 
and Folger, 1993]. Some of the relevant aspects of fuzzy set theory are presented in 
Appendix K. 

In this chapter a parametric study is carried out to study the influence of various 
parameters on the array size determination. The parameters which play dominating 
role, are taken as fuzzy. The possibility distribution of array size is obtained for different 
months. If an array is designed on the basis of worst expected conditions, it will be 
a reliable design. But the cost of the airconditioning system will be exorbitant. On 
the other hand if it is designed on the basis of best expected conditions, the cost will 
be much less but reliability will be poor. Therefore, it is prudent to design at some 
intermediate operating conditions which provide an optimal solution. A method to 
select an optimum array size obtained from the fuzzy solution is discussed. 

3.2 Parametric Study 

The array size is dependent on the cooling load and power generated by photovoltaic 
array. Cooling load and array output are dependent on a number of factors. Some of 
these factors play dominating role, while other influence only slightly. A parametric 
study is done to find out the effect of various parameters, so as to enable us to find out 
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critical parameters. These parameters will then be taken as fuzzy parameters. The 
output for different values of parameters can be found by changing the input variables 
in the input file of cooling load estimation. 

(a) Effect of Maximum Temperature: From Figure (3.1), it is clear that the maximum 

temperature of the day has much influence on cooling load. A 10°C tempera- 
ture increase, causes cooling load to increase by a factor of 1.7. However, from 
Figure (3.2), it is seen that array output increases insignificantly with increasing 
maximum temperature. 

(b) Effect of Minimum Temperature: Figure (3.3), shows that like maximum temper- 

ature of the day, minimum temperature of the day also influences cooling load. 
Increase in the minimum temperature of the day increases cooling load. However, 
Figure (3.4), exhibits that increase in array output is not significant. This aspect 
should be kept in mind while adopting the use of photovoltaic airconditioning 
system. 

(c) Effect of Humidity: Figure (3.5) shows an increase in cooling load with increasing 

humidity. Humidity associated with ventilation and infiltration air increases the 
cooling load only. But humidity does not have any influence on array output. 

(d) Effect of clouds: Figure (3.6) shows that increasing cloud cover decreases cooling 

load only slightly. This is because firstly with clouds only beam radiation reduces, 
diffused radiation being slightly increased and secondly since the most of the 
portion of the room is already covered, radiation in itself accounts for small 
fraction of cooling load. However, from Figure (3.7), it is seen that increasing 
cloud cover decreases array output drastically. 
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(e) Effect of Number of Occupants: It is seen from Figure (3.8), that as the number 

of occupants increases, cooling load also increases. 

(f) Effect of Sunrise and Sunset time: The sunrise and sunset times on the first and 

the last day of month are different. However, due to change in sunrise and sunset 
timings, cooling load and array output do not change much. Therefore, taking 
average values of sunrise and sunset timings in a month is appropriate. 

(g) Effect of Door Openings: It is observed that cooling load during normal and 
heavy door openings are not much different. Prom the above discussion it is clear 
that critical parameters are: maximum and minimum temperature of the day, 
humidity, cloud cover and number of occupants. 


3.3 Fuzzy Set Modeling of Various Parameters 


Maximum and minimum temperatures of the day, humidity, number of occupants in 
the room and fractional cloud cover are five parameters which have been considered 
fuzzy. These parameters are characterized by their membership functions. There are 
various ways to construct a membership function [Vallipan and Pham, 1995], however, 
for simplicity only linear triangular membership function is considered here. Triangular 
membership function is given by. 


IJ.{x) = 


0 

X < 1', 

x-l' 

m-V 

1' < X <m, 

h* —X 

h^—m 

m < X < h', 

0 

X < h' 


(3.1) 


where 
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(3.2) 


and 


[ 0 m < 2{m — 1) 

h' — m + 2{h — m). (3.3) 


Here /, m and h are experts estimates of a parameter for low, most likely, and high values 
respectively. The value of // is taken to be 0.5 at x = I and h and 1 at rr = m. These 
fuzzy parameters lead to the fuzzy arithmetic computations at each a-cut. Design 
parameters can be obtained then as fuzzy numbers. 

For each parameters, different membership functions are constructed for different 
months. Average values of each month are considered most likely values, while 0.9 
times average value is considered low and 1.1 times average value is considered high. 


3.4 Results with Fuzzy Parameters 

The programme was run for estimation of PV array size for airconditioning a room in 
Kanpur. The particulars of room are given in Appendix G. COP of the airconditioning 
system is taken as 3.84 based on section 7.1.3 and Chapter 7 of reference [Prasad, 
1993]. PV cell data is given in Appendix H. Figures (3.9-15) show the possibility 
distributions of cooling load for different months (from April to October). There is 
a great variation in the cooling loads of various months. Figures (3.16-22) show the 
possibility distribution of array wattage required for different months. It is to be noted 
that there is not a direct correlation between cooling load and wattage of array, as the 
cloud cover is having much influence on the array output. 

The question arises that given possibility distributions of array wattage for differ- 
ent months, which wattage of array should be selected. This is a subjective decision. 
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There can be various methods for it. One method is to give linguistic expressions to 
membership grades. For example, if the array is selected on the basis of upper limit 
of 0.5 membership grade, the design performance can be called excellent, because it 
is based on the worst weather conditions. The design based on the upper limit of 0.6 
membership grade is very good and those based on the upper limit of 0.7, 0.8 and 
0.9 will be called good, satisfactory and poor respectively. The design based on the 
membership grade of 1, will be called very poor. Once a particular array is selected, it 
will fall in the different categories of design in different months. To assess the overall 
effectiveness of the design, grade points are allotted to design categories as follows; 


Design category 

Grade point 

Excellent 

10 

Very good 

8 

Good 

6 

Satisfactory 

4 

Poor 

2 

Very poor 

0 


Then the average grade point for all months can be found and is converted back 
into linguistic expression. 

The procedure is similar to the evaluation of overall performance of a student in all 
the courses, he has taken. He is awarded a letter grade in each course, and each grade 
has corresponding points. To evaluate his overall performance, cumulative performance 
index is obtained by adding all the points, he has obtained in his various courses and 
dividing it by the number of courses, he has taken. Similar to that we obtain here, 
average grade point for the photovoltaic array, over the various months for the array. 
Examples: 

Suppose an array of 4500 W is selected. This design will be satisfactory in the 
month of August and in the rest of months, from April to October, it will be excellent. 




(3.4) 


The average grade point is given by, 


av. g. p. = 


10 + 10 + 10 + 10 + 4 + 10 + 10 
7 


9.14 


which is between excellent and very good design. 

Suppose an array of 3000 W is selected. This design is excellent in April, May, 
June and July, very poor in the month of August, poor in the month of September and 
excellent in the month of October. Therefore, 


av. g. p. = 


10 + 10 + 10 + 10 + 0 + 2 + 10 

7 


= 7.43 


(3.5) 


This shows that this design is between good and very good, more inclined towards 
very good design . The cost of 4500 W array is about 1.5 times that of 3000 W array. 
Therefore a designer may decide to choose 3000 W array, as the quality of design is 
not drastically sacrificed. 

As a last example, suppose a 2300 W array is selected which gives satisfactory 
performance in the month of May (month of maximum heat). This design is excellent 
in the month of April, good in the month of June, satisfactory in the month of July, very 
poor in the months of August and September and excellent in the month of October. 
Then, 


10 + 4 + 6 + 4 + 0 + 0 + 10 
av. g. p. = 


= 4.86 


(3.6) 


which is between satisfactory and good design. This design may be opted where cost 
is the primary consideration. 


3.5 Summary 

A method to estimate array size for airconditioning the space, using fuzzy set theory, 
is presented and computer code has been developed for this purpose. An example has 
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been taken for estimating array sizes for a typical room at IIT Kanpur. It is highlighted 
that use of fuzzy set theory in the array size determination may be much helpful in 
selecting a optimum array size. 
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Chapter 4 


OPTIMUM LOAD MATCHING 
OF A BATTERY CONNECTED 
PHOTOVOLTAIC ARRAY 

4.1 Introduction 

In a photovoltaic system it is desirable to extract the maximum amount of energy 
output from the array connected to the battery having proper matching. This matching 
is usually done by the suitable choice of number of series modules and number of parallel 
strings of the photovoltaic array. In addition to this battery parameters may also be 
decided suitably in order to have maximum energy transfer. 

In the past researchers [Khouzam, 1990; Appelbaum, 1987], have concentrated on 
optimum array for photovoltaic power. In this chapter, a method to select optimum 
layout for PV cells and battery parameters is presented. The importance of proper 
load matching for airconditioning and refrigeration is highlighted. 

In the present work, the photovoltaic power generated by the array at any instant is 
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considered to be proportional to the solar radiation falling on it, unlike the assumption 
of sinusoidal variation of the photovoltaic current with solar time, made by earlier 
researchers [Khouzam, 1990; Appelbaum, 1987]. A comparison with the experimental 
results validates our model. Using this model a parametric study is carried out. 

Alternatively, to extract the maximum energy output from the PV array, a max- 
imum power point tracker (MPPT) can also be employed in between PV array and 
the load. The maximum power point tracker is an electronic circuit which can be used 
to operate the load at the maximum power point of the array. However, MPPT is a 
costly and use of an additional circuit needs frequent maintenance. Therefore, attempt 
should be made to avoid its use by proper load matching as discussed here. 

4.2 Load Characteristics 

The battery is directly connected to PV array. The load equation for the battery is 
given as 

Vl = I.R, + V, (4.1) 

where Vi is the load voltage, I is the load current, 14 is the induced emf of the battery 
and Re. is the internal resistance of the battery. If the small potential drop from array 
to battery is neglected, then 


Vl = V (4.2) 

where V is the output voltage of the PV array. Substituting the expression for V from 
equation (4.1) into the equation (2.1), one gets 


Ve = AVt. ln{ilph -1 + Ir)/Ir} " ^Rs + Re) 


(4.3) 



The power stored in the battery is given by 

Pb = I.Ve = AVT.Ihx{{Iph -I - /r)//r} “ P{Rs + Re) (4.4) 

while the power available at the array is given by 

Pa = PV = AVT.I\n.{{Iph — I + 7^)//^} — PRg (4.5) 

The total daily energy which can be stored in the battery, assuming that the battery 
is connected to the array during the entire period, is then calculated from : 

rtss 

Eb — Ppdt ( 4 :- 6 ) 

tsr 


4.3 Load Matching Factor 

Each array possesses a maximum power line. If the array operates at maximum power 
point voltage at all the time, the energy output from the array will be maximum. De- 
noting the maximum possible energy storage in the battery by Epmax, a load matching 
factor is defined as 

M = Ep/EBmax (4.7) 

The closer is the value of M to 1, the better is PV array matching. 

The procedure to calculate Epmax is described here briefly. The maximum power 
is found by differentiating power (P) in equation (2.9) with respect to load current I 
and equating the same to zero. This gives the following equation: 

Iph = Imp + 7r {exp[(2/mpi?s/AVT) + Imp! {Iph “ Imp + 4)] “ 1} (4.8) 

This equation is solved numerically to find Imp- 

The maximum power point voltage, Vmp can be found from equation (2.1) by re- 
placing I by Imp- The maximum power, Pmp is then given by: 

Pmp = Vmp-Imp (4-9) 
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The maximum electrical energy which can be generated from the solar energy over one 
day period is given by the integral: 

ptss 

^Bmax / ( 4 . 10 ) 

'Jtsr 

which is a characteristics of the array. Here, t^r and t,, are the sun-rise and sun-set 
times respectively. 

4.4 Selection of Battery Voltage and Capacity 

It is observed that most of the refrigerators require fractional horse power motors. A 
permanent magnet D.C. motor powered by 12 volt battery source may be employed as 
the D.C. motor is inherently efficient for part load operation. For airconditioners and 
heaters requiring high power motors either brushless D.C. motor or induction motor 
along with optimal flux control circuit (to make it part load efficient) may be employed. 
For this purpose a 24-volt battery along with inverter may be used. 

The battery capacity is decided from the consideration of maximum reserve energy 
which we want to store and maximum depth of discharge. Usually a battery capacity 
of twice the total solar energy of the day will be sufficient. However, where the weather 
conditions are uncertain a higher battery capacity should be used. 

4.5 Problem Formulation 

The objective is to maximize the load matching factor and consequently Eb, since 
Esmax is fixed for a given array and insolation and temperature profile. This can be 
stated mathematically as follows: 

Maximize: 

EB = f{Ns,Np,Re) (4.11) 
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Subject to: 


Ns.Np < Na (4.12) 

where Na is the maximum number of cells in the array. 

The solution of the above optimization problem yields the optimum storable battery 
energy, from which the optimum load matching factor, Mopt and the corresponding 
series cells, parallel strings and internal resistance can be calculated. 

4.6 Solution Procedure 

The optimization problem is solved by the sequential quadratic programming method 
[Pike, 1986]. The E04 UCF routine of NAG library is used for this purpose. In this 
routine the design variables Ns, Np and Re are treated as real numbers. After the final 
solution is obtained, the nearest integer value of Ns and Np is chosen. 

In equation (4.3) I and in equation (4.8) Imp is found using interval reducing meth- 
ods [Pike, 1986]. The interval reducing methods guarantees the convergence of the 
solution. The overall computation time for the problem is 1-2 seconds of user time on 
HP-computers. 

4.7 Results and Discussion 

The PV generator data are obtained from literature [Appelbaum, 1986]. The two array 
sizes were selected - one array consisting of 400 cells is suitable for 100 W refrigerator, 
the other array consisting of 7040 cells is suitable for airconditioning of room. The 
results obtained from the simulation runs are discussed in the following subsections. 
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4.7.1 Optimized Combination of Cells and Internal Resistance 

Xabl6 4.1 shows tho optimizod p&r&nietors for tho array and battory. It is obsorved 
that in all the cases internal resistance has to be kept zero. Physically zero internal 
resistance is not possible, but attempt should be to keep it at minimum. The optimum 
energy stored in the battery was 0.98 times maximum storable energy EBmax for bigger 
array and 0.97 times maximum storable energy Esmax for smaller array. 

Here the optimization was done for maximizing the stored energy. If the optimiza- 
tion is done to maximize available energy the different solutions are obtained. For 
example, for 24-V battery the smaller array will give maximum available power for Np 
= 6, Ns = 66 and = 0.5 fl; however, then a fraction of the energy will be wasted 
as the heat loss and useful energy will be only 0.92 times Esmax- Since our concern is 
to maximize stored energy, the optimization has to be done for maximizing the stored 
energy and not the available power. 

Figure (4.1) shows the effect of proper array layout on the load matching factor. 
The curve has been drawn for 12-V battery and 400 cell array combination. On the 
horizontal axis, number of parallel strings are shown. The number of cells in series 
may be found by taking the nearest integer value of 400/Np. The Figure clearly shows 
the proper cell arrangement is the crucial thing to decide. The load matching factor 
is maximum at Np = 12, and more the deviation from this value, we have, the lesser 
load matching factor will be obtained. 

4.7.2 Effect of Internal Resistance 

To study the effect of internal resistance on the load matching factor, first a system 
consisting of 12- volt battery and 400 cell array is considered. In Figure (4.2), the lower 
curve is for the load matching factor based on stored power, while the upper curve is 
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Table 4.1: Optimized Parameters for the Array 



Smaller array 
(400 cells) 

Bigger array 
(7040 cells) 


Np Ns 

re 

Np Ns 

re 

For 12-V 
Battery 

12 33 

0 

227 31 

0 

For 24-V 
Battery 

6 66 

0 

117 60 

0 


for available power. The figure shows clearly that internal resistance has the adverse 
effect on the load matching factor based on the stored power. However, observing the 
load matching factor based on the available array output, it is seen that it attains 
maxima at some value of internal resistance. Figure (4.3) is for bigger array (7040 cell 
array). Here, also the same trend is observed, with different pattern. 

In practice, there is always some resistance in the battery. The question arises 
whether the array designed for the battery with zero internal resistance will be opti- 
mum for the battery with the resistance. From Table 4.2, we see that answer is No . 
For example, for 1.2Q internal resistance battery, the optimum array arrangement is 
Np = 9, iV^ = 44 and available array energy is 0.98 Esmax which shows if we design the 
array suitable, the maximum power point tracker (MPPT) is not needed. However, 
because of i'^R losses, the stored energy is less. On the other hand, if we take Np = 
12, Ns = 33, (which is optimum for 0 internal resistance battery) for 1.20, internal 
resistance battery we obtain lesser available and stored energy. This shows that with 
the change of internal resistance optimum combination also gets changed. 
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Table 4.2: Arrangement of smaller array(400 cells) for 1,20 internal resiste npR battery 
Re Stored energy Available energy Maximum energy 

1.2 n 9 44 636.16 kWh 922.71 kWh 928.6 kWh 

1.2 0, 12 33 474.38 kWh 633.61 kWh 928.6 kWh 

4.7.3 Sensitivity Analysis with Changing Battery Voltage 

Figures (4.4-5) show the influence of battery voltage on the load matching factors for 
the arrays designed for 12 volt and zero internal resistance battery. It is seen that 
load matching factor drops first gradually, then at the high rate as the battery voltage 
deviates from 12 volt. The attempt should be to keep the variation of the battery 
voltage to the minimum possible value. 

4.7.4 Load Matching Factor for Different Months 

The load matching factors for different months for the smaller and bigger arrays were 
found. It was s(Km that the load matching factor is almost same for all months. This 
shows that the array which has been matched in one month will continue to operate 
at almost maxinmm load matching factor in other months as well. 

4.7.5 Summary 

In this chapter, an optimization code has been developed and simulation study is done. 
Based on the study, the following conclusions may be drawn. 

1. There is no need of a maximum power point tracker (MPPT) if the array is 
properly matched with the battery. Because the value of load matching factor in 
the range of applications turn out to be around 0.97. 


2. The internal resistance of the battery has to be kept at minimum. 





0.95 



Fig. 4, 5 Variatio 
(np = 22' 



3. In practice the internal resistance cannot be kept zero. Therefore, the array 
should be matched foi Np and keeping the specified internal resistance value 
constant. 

4. If an array is matched in one month, it will give good results in other months 
too. 
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Chapter 5 


ARRANGEMENT OF PV ARRAY 
AND ENVIRONMENTAL 
BENEFITS 

5.1 Introduction ^ 

The design of a photovoltaic array to meet the cooling load requirement of a room 
is very vital from the vi(nv point of economy. Once an array is selected, it has to be 
physically arrangcal in order to g(d maximum output from it. Usually, an array is put 
on the top of a roof. Tlum, the radiation falling on the roof is fruitfully converted 
into electrical energy. In this chapter, a complex search algorithm is used to reduce the 
array size by its optimum arrangement and tilt angle to get maximum power for known 
photovoltaic characteristics and cooling load. This is done to reduce the array size to 
as low value ;is i)()ssible. The a<ioption of photovoltaic array eliminates the radiation 
energy du('. to covertal roof, rendering reduced cooling load and hence reduced power 
for airconditioning, romputations were carried out for replacing electrical energy from 
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the conventional fuels by photovoltaic cells and subsequent reduction in global warming 
gases like C02- Results are presented in the graphical form. Thus photovoltaic arrav 
plays following roles: 

• The reduced cooling load and hence reduced power for airconditioning. 

• Elimination of production of global warming gases and thus reduced global warm- 
ing effect. 

• Conservation of depleting fossil fuels. 


5.2 Reduction in the Cooling Load due to Coverage 
of Roof 

Heat transfer through walls and roofs is a major component of cooling load. The solar 
radiation may render the roof temperature as high as 355 K in an environment at 315 
K. Hence, installing the photovoltaic cells over roof provides a cover to the roof. The 
array size is found to be much more than the size of the roof. By mounting the array 
on the terrace, therefore, cooling load requirement is reduced because of the coverage. 
Figure (5.1) shows the array mounted on the terrace. A code has been developed which 
computes cooling load on hourly basis as given in Section 2.4. 

The hourly cooling load helps determine a suitable size of the airconditioning sys- 
tem, and components can be designed accordingly. Integrating cooling load with re- 
spect to time using Simpson’s formula, the energy requirement for 24 hour period is 
estimated, which is used for the estimation of PV array size. Using this code, one can 
find the cooling load requirement of an uncovered room and an array-covered room. It 
will be seen that a large portion of cooling load requirement is met just by covering 
the roof. 
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Fig. 5,1 PV array covered roof. 




5.3 Optimum Arrangement of Array 

From Section 2.2 and Section 2.3, it is clear that power output from the array is depen- 
dent on the slope /? of the array, which is the angle made by the plane surface with the 
horizontal and the surface azimuth angle 7 which is the angle made in the horizontal 
plane between the line due south and the projection of the normal to the surface on 
the horizontal plane. Two optimization routines (successive quadratic programming 
and complex search method) [Deb, 1995] were used for finding the optimum ^ and 7. 
The successive quadratic programming takes much time, because it has to calculate 
gradients by finite difference technique, while in complex search method gradient infor- 
mation is not needed. The algorithm of complex search method is given in Appendix 
L. 

The optimization problem is; 

Maximize jE'b(/?, 7) (5-1) 

0 < /3 < 90, -90 < 7 < 90 (5.2) 

where Eb{/3, 7) is the energy stored in the battery as given in the equation (4.7) which 
is dependent on the variables /? and 7, besides other constant parameters. Sometimes 
the slope of the array may have to be fixed from the consideration that in rainy season, 
water falling on it should flow smoothly. In that case only 7 is optimized, the value of 
/? being a fixed value. 

5.4 Results and Discussion 

For parametric study a 7040 cell array is used which is suitable for airconditioning of 
a room details of which are given in Appendix-G. The PV generator data are given in 
Appendix-H. 
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5.4.1 Energy Saving due to Roof Coverage 

As described in Section 2.4, a code for computing cooling load on hourly basis has 
been developed. This code was run for the room data given in Appendix-G. Figure 
(5.2) shows the cooling load variation with time for the month of May for a 3-side 
shaded room, in which the south side is exposed to the sun. Results for the covered 
and uncovered roof are shown in the same figure. It is seen that most of the time, 
the cooling load is different for the covered and uncovered roof rooms. Actually, the 
radiation falls only at day time, but due to time lag, its effect is experienced upto 
late night. The total cooling load for the uncovered roof is 72.22 kWh and for the 
covered roof it is 60.48 kWh. This shows that we obtain 16% saving in the cooling load 
by covering the roof. Figure (5.3) shows the variation of cooling load with time for 
the month of August. Figures (5.4-5) show electrical power needed for airconditioning 
respectively in the month of May and August. Figures (5.6-7) show the power given by 
the PV array as well as saving in electrical energy due to roof coverage. The effective 
power is the sum of two energies and is also shown in Figures (5.6-7). 

Variation of cooling load with time for the same room conditions except that all 
the walls are exposed to the sun is exhibited in Figure(5.8). In this case also, the 
saving in the cooling load due to roof coverage is observed clearly. The total cooling 
load in a day is 83.74 kWh for uncovered roof and 72.01 kWh for covered roof. This 
renders about 14% saving in the cooling load with the covered roof. Figure (5.9) shows 
the cooling load variation for the month of August. Figures (5.10-11) show electrical 
power needed for May and August. Figures (5.12-13) depict the power obtained due 
to photovoltaic array. 

Thus, a significant amount of radiation energy is saved by the array with its cov- 
erage. It can, therefore, be concluded that the large space required by the array is 
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Fig, 5. 6 Effective power output from PV array v/s Time. 
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Fig -5. TO Electric power needed for Airconditioning v/s Time. 










Fig -5.13 Effective Power output from PV array v/s Time 




not merely a drawback of PV array. It has its positive effect too. In fact in an unsur- 
rounded room (roof and walls exposed to the sun), the roofs as well as walls can be 
covered by the large size of PV array required for airconditioning that room. In the 
case of example room, the cooling load for a bare room is 83.74 kWh for a day, whilest 
it is 60.19 kWh in a day, for a fully covered room. Thus about 28.12% cooling load is 
reduced in this case. 

5.4.2 Optimum Arrangement of the Array 

For Kanpur city, the optimum slope /9 comes out to be zero. The surface azimuth angle 
7 can be any value. However, array is usually kept at some angle in order to avoid free 
rain-fall. It is found that if array is kept at 5° angle, the total energy obtained is only 
0.42% less than that obtained from horizontal array, whilest if it is kept at 15° angle, 
energy is 3.4% less. It is found that if the array is kept at some slope, 7 should be 90°. 

If the latitude of the place is 50°, then P = 18.7° and 7 = 0 ° combination gives 
maximum power, whilest if the latitude is 80°, then = 46° and 7 = 0 ° combination 
gives maximum power. Thus it is seen that latitude of a place has much influence on 
the optimum arrangement of the array. 

5.4.3 Reduction in the Atmospheric Pollution 

When the electricity is generated by conventional fuels, like coal, oil or natural gas, 
CO 2 is produced. Figures (5.14- 15) show how much CO 2 will be reduced if photo- 
voltaic array replaces these sources. The data about CO 2 production is taken from 
[HR Bulletin, 1994]. These values are 1 . 12 , 0.94 and 0.57 kgs of C 02 /kWh for coal, oil 
and natural gas, respectively. 

In addition to this, there are transmission and distribution losses which are taken 
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Fig* 5. 18 Kg of Coal-ash saved Vs time 




to be 25% on all India basis as given in [UPSEB, 90-91]. Therefore emission of CO 2 is 
further increased by 25% due to these losses. 

It is seen that for airconditioning of a room size as given in Appendix G, CO 2 at the 
rate of as much as 3.75 kg/s at a peak time can be saved as shown in Figures (5.16-17). 
If photovoltaic energy is used in place of electricity, then ash which is produced when 
coal is burnt to generate electricity is eliminated. This solves the ash disposal problem 
in addition to rendering cleaner environment. The amount of ash saved is shown in 
Figure (5.18). 

5.5 Summary 

In this chapter energy saved due to roof coverage is discussed. The same is concluded 
as under: 

• A complex search method code is used to find out the optimum arrangement of 
the array. 

• It is seen that this arrangement is dependent on the place. 

• It is shown that the array not only supplies the energy, but also reduces energy 
needed by providing shading. 

• Finally, the reduction in the CO 2 in the atmosphere is highlighted. 
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Chapter 6 


CONCLUSIONS AND 
SUGGESTIONS 

6.1 Conclusions 

6.1.1 Estimation of Photovoltaic Array Size for Aircondition- 
ing, Heating, and Refrigeration 

A deterministic method to estimate array size for airconditioning heating and refrig- 
eration system is presented and computer code has been developed for this purpose. 
For each of these cases, examples have been taken for estimating array sizes. For 
airconditioning and heating both vapour-absorption and vapour- compression systems 
were considered. It is observed that array size requirement is much more in the case 
of vapour-absorption system than in the case of vapour-compression system for both 
airconditioning and heating. This is due to the low COP of the vapour-absorption 
system. 

For the case of a refrigerator, 100 W, 45 L capacity refrigerator can be constructed 
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which can successfully operate at Kanpur and Bombay, the two representative cities of 
India where environmental conditions differ much from each other. 

6.1.2 Application of Fuzzy Set Theory for Array Size Esti- 
mation 

Fuzzy set theory is employed with a view to get more realistic estimate of array size 
for airconditioning the space. A computer code has been developed for this purpose 
and an office room at IIT Kanpur has been selected for estimating array size for the 
same. The use of fuzzy set theory in the determination of array size and selecting an 
optimum size of array is highlighted . 

6.1.3 Optimum Matching of the Array with Battery 

Direct matching of photovoltaic array with battery eliminating the need for maximum 
power tracker is suggested. An optimization code has been developed and simulation 
study is done. Based on the study, the following conclusions may be drawn: 

• There is no need of a maximum power point tracker (MPPT) if the array is 
properly matched with the battery. Because the value of load matching factor in 
the range of applications turn out to be around 0.97. 

• The internal resistance of the battery has to be kept at minimum. 

• In practice the internal resistance cannot be kept zero. Therefore, the array 
should be matched for Ns and Np keeping the specified internal resistance value 

constant. 

• If an array is matched in one month, it will give good results in other months 
too. 
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6.1.4 Saving of Pollution by using Photovoltaic Array 

• Energy saved due to roof coverage is discussed. 

• A complex search method code is used to find out the optimum arrangement of 
the array. 

• It is seen that this arrangement is dependent on the place. 

• It is shown that the array not only supplies the energy, but also reduced energy 
needed by providing shading. 

• Finally, the reduction in CO2 in the atmosphere is highlighted. 

A complex search method code is used to find out the optimum arrangement of the 
array. It is seen that this arrangement is dependent on the place. It is shown that 
the array not only supplies the energy, but also reduced energy needed by providing 
shading. 

6.2 Suggestions 

The present investigation has given a lot of scope for further work. The same is 
presented below: 

6.2.1 Development of Photovoltaic Airconditioner without 
Battery: 

This is suitable for airconditioning an office area during the day time. For this a 
maximum airconditioning tracker (MAT) should be developed which operates the pho- 
tovoltaic array at an operating point (V*, T) corresponding to which maximum aircon- 
ditioning will be achieved. MAT will operate the PV array at (V*, F) which results in 
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maximum rpm of the compressor motor. Maximum rpm of compressor motor results in 
maximum airconditioning as H.P. of the compressor is proportional to its rpm. MPPT 
is provided in between the PV array and the airconditioner motor. 

This will be ideal for airconditioning offices. Provision will also be made to run the 
airconditioner from the electric supply available in the office. Therefore, if the need 
be, airconditioner can also be operated at night. 

6.2.2 Development of Photovoltaic Airconditioner with Bat- 
tery: 

This airconditioner is ideal for a place where there is no power supply or where air- 
conditioning is required for all the 24 hours without interruption caused by power 
failures. 

In this version, a battery is provided after the PV array. Energy from photovoltaic 
array is stored in the battery and the airconditioner is run from the battery. 

6.2.3 Development of Photovoltaic lOOW Refrigerator: 

This refrigerator will be ideally suitable for storing medicines in dispensary in villages 
where there is no power supply. 

6.2.4 Development of Electrically Controlled Compressor: 

An electrically controlled multi-cylinder compressor to eliminate the problem associated 
with its starting (High Starting Current) should be developed, which can be used in 
photovoltaic airconditioner/refrigerator of high tonnage capacity to enhance energy 
conservation by providing load matching between load (number of loaded cylinders 
can be varied electrically) and photovoltaic array. 
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A multi-cylinder compressor with a capability of operating with one, two or three 
cylinders unloaded electrically can be developed. In the presently available multi- 
cylinder compressors, the unloading is achieved by opening the suction valves of the 
unloaded cylinders by some mechanical mechanisms driven by the suction pressure. In 
the electrically controlled multi-cylinder compressor, unloading should be achieved by 
using electrically actuated devices which will permit the compressor stage unloading 
operation to be carried out in response to electrical control signals. This method will 
not require any suction pressure. A simple controller will be developed to vary the 
number of stages as the solar insolation varies automatically unloading more stages as 
the insolation decreases. This will enhance energy conservation be improving matching 
between load and PV array. 

With this development, it will be possible to start any airconditioner /refrigerator 
with all cylinders unloaded thereby having low starting current for the compressor 
motor. This will eliminate the electrical dip associated with starting of any aircondi- 
tioner / refrigerator. 
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Appendix A 


Values of the Constants A, B and C 
used for Predicting Hourly Solar 
Radiation on Clear Days 


Sukhatme, 1984 



A(W/m) 

B 

C 

January 21 

1228 

0.142 

0.058 

February 21 

1213 

0.144 

0.060 

March 21 

1185 

0.156 

0.071 

April 21 ■ 

1134 

0.180 

0.097 

May 21 

1103 

0.196 

0.121 

June 21 

1087 

0.205 

0.134 

July 21 

1084 

0.207 

0.136 

August 21 

1106 

0.201 

0.122 

September 21 

1150 

0.177 

0.092 

October 21 

1191 

0.160 

0.073 

November 21 

1219 

0.149 

0.063 

December 21 

1232 

0.142 

0.057 
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Appendix B 


Transmissivity and Absorptivity of 
Ordinary Glass 


[Arora, 1993] 


Angle of Incidence Transmissivity Absorptivity 
(in degrees) 

0 

0.87 

O.Of) 

20 

0.87 

0.05 

40 

0.86 

0.06 

50 

0.84 

0.06 

60 

0.79 

0.06 

70 

0.67 

0.06 

80 

0.42 

0.06 

90 

0.00 

0.00 


133 




Appendix C 


Infiltration load 

[Prasad, 1993] 
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Room A umber oi air rhangt>; 
'■olume per 24 h 

(m^) Average load 

~5 ’22.G ’ 

10 20.5 

15 19.25 

20 18.25 

25 17.37 

30 16.60 

35 15.91 

40 15.29 

45 14.74 

50 14.25 

55 13.80 

60 13.40 

65 13.03 

70 12.70 

75 12.40 

80 12.12 

85 11.87 

90 11.67 

95 11.43 


por 2 } h 
heavy hiufl 

I, i ; 1 

32,2(5 

30.14 

28.7(5 

27.27 

2(5.04 

25.10 

24.45 

24.08 

23.99 


dumber of air Room volume m 
changes per 

. h long usage long usage 


12.65 

55 

11.26 

60 

10.69 

70 

9.92 

80 

8.92 

90 

8.26 

100 

7.22 

120 

6.44 

140 

6.13 

150 

5.85 

160 

5.37 

180 

4.99 

200 

4.42 

240 

4.30 

250 

4.01 

280 

3.86 

300 

3.55 

350 

3.38 

400 

3.19 

450 

3.08 

500 

2.94 

600 

2.86 

700 

2.82 

800 

2.80 

1000 

2.80 

1500 

2.80 

2000 

2.75 

2500 

2.73 

3000 


Appendix D 


Ventilation air requirement 

[Prasad, 1993] 
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Application 

Smoking 

Ventilation 

(m^/s per person)x 

10-3 

Recommended 

Minimum 

Ventilation 
m^/s/m^ of 
floor area 
Minimum x 
1000+ 

Apartment, average... 

Some... 

9.5 

4.7 

1.0 

Apartment, Deluxe... 

Some... 

14.2 

12.0 

1.676 

Baking space... 

Occasional 

4.7 

3.5 


Barber shops... 

Considerable 

7.1 

4.7 


Beauty parlours... 

Occasional 

4.7 

3.5 


Brokers’board rooms. 

Very heavy 

24.0 

9.5 


Cocktail bars... 

Heavy 

19.0 


12.0 

Corridors (Supply or 
exhaust) 



_ 

— 

1.27 

Department stores 

None 

3.5 

2.4 

0.254 

Directors’ room 

Extreme 

24 


14.2 

Drug stores+ 

Considerable 

4.7 

3.5 


Factories** 

None 

4.7 

3.5 

0.508 

Five and ten cent 

stores 

None 

3.5 

2.4 


Funeral parlours 

None 

4.7 

3.5 

5.08 

Garages++ 

- 

— 

— 

Hospitals operating 
rooms 

None 



10.16 

,, private rooms 

None 

14.2 

12.0 

1.676 

,, wards 

None 

9.5 

4.7 


Hotel rooms 

Heavy 

14.2 

12.0 

1.676 

Kitchens, Restaurant 

- 

- 

— 

20.32 

,, , Residence 

- 

- 

- 

10.16 

Laboratories+ 

Some 

9.5 

7.1 

6.35 

Meeting rooms 

Very heavy 

24 

14.2 

Offices, General 

Some 

7.1 

4.7 

1.27 

Private 

None 

12.0 

7.1 

Private 

Considerable 

14.2 

12 

1.27 

Restaurant cafeteria+ 

5 

,5.7 

4.7 


Dining Room+ 

57 

7.1 

0.7 


School Rooms 

None 

— 



Shop, Retail 

None 

4.7 

3.6 


Theatre** 

None 

3.6 

2.4 ■ 


Theatre 

Some 

7.1 

4.7 

10.16 

Toilets** (exhaust) 

- 

- 

— 
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* When minimum is used take the larger of two. 

See local codes which may govern. 

+ May be governed by exhaust. 

++ May be governed by special sources of contamination of local codes. 
All outside air recommended to overcome explosion hazard of anesthetics. 
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Appendix E 


Heat Liberated due to Occupancy 


[Arora, 1993] 



Metabolic 



Heat Liberated, W 



Room Dry Bulb °C 

Activity 

Rate 

W 

20 

22 

24 


26 



S 

L 

S L 

S 

L 

S 

L 

Seated at rest 

115 

90 

25 

80 35 

75 

40 

65 

50 

Office work 

140 

100 

40 

90 50 

80 

40 

70 

70 

Standing 

150 

105 

45 

95 55 

82 

68 

72 

78 

Eating in 









restaurant 

160 

110 

50 

100 60 

85 

75 

75 

85 

Light work in 









factory 

235 

130 

105 

115 120 

100 

135 

80 

155 

Dancing 

265 

140 

125 

125 140 

105 

160 

90 

175 


S= Sensible Heat 
L= Latent Heat 
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Appendix F 


Heat of Respiration of Products in 
J /kg per 24 hours 


[Arora, 1993] 



Product 


Storage Temperature 



0°C 

4.4°C 

15.6‘’C 

Apples 

312-1560 

625-2810 

2390-8215 

Bananas 


— 


Cabbage 

1248 

1770 

4265 

Carrots 

2183 

3640 

8420 

Cauliflower 

- 

4680 

10500 

Cherries 

1352-1871 

- 

11400-13725 

Cucumbers 

- 

- 

2290-6860 

Grape fruit 

416-1040 

730-1350 

2290-4160 

Grapes, American 

624 

1250 

3640 

Grapes, European 

312-416 

- 

2290-2705 

Lemons 

520-936 

625-1975 

2390-5200 

Melons 

1350 

2080 

8840 

Mushrooms 

6446 

- 

- 

Onions 

728-1144 

830 

2495 

Oranges 

416-1040 

1350-166 

3850-5405 

Peaches 

936-1456 

1455-2080 

7590-9670 

Pears 

728-936 

- 

9150-13725 

Peas 

8526-8733 

13520-16635 

40860-46265 

Plums 

416-728 

935-1560 

2495-2910 

Potatoes, immature 

- 

2705 

3015-7070 

Potatoes, mature 

- 

1350-1870 

1560-2705 

Strawberries 

2807-3950 

3745-7070 

16220-21105 

Tomatoes, green 

625 

1145 

6445 

Tomatoes, ripe 

1040 

1350 

5820 

Turnips 

1975 

2290 

5510 
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Appendix G 


Specifications of the room selected 
for study 


Room size = 4.85 x 4.30 x 4.20 
Thickness of the wall = 0.24m 
Thickness of the roof = 0.12m 
Number of windows = 3 (South side) 

Area of the window = 7.2 

Only south side exposed to the sun. 

Number of occupants = 4 
Number of fans = 1 
Number of lights = 2 
Ratings: 

Fan = 60 Watts 
Light = 40 W 

The walls of the room have layers of plaster of 15 mm on both sides. The conductivity 
of plaster is 0.72 W/m°C and that for brick structure 1.30 W/m°C. The roof is made 
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of concrete with 15 mm mortar layers on both sides. The conductivities of mortar and 
concrete are 0.72 W/m°C and 1.73 W/rrfC respectively. The absorptivity of whole 
structure is taken to be 0.3. Time lag and decrement factor for the room is taken from 
reference [Arora, 1993]. 
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Appendix H 


Specifications of Photovoltaic Array 

[ Appelbaum, 1986] The PV array data per cell at 100 % sun and 25°C used for the 

simulation are as follows: 

Iq = 0.8 Amp. 

Ir = 0.0005 Amp. 

Vt = 0.731 Volt 
R, = 0.05 Ohm 
a = 5.70 X 10-"^ 
b = 4400 
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Appendix I 


Particulars of the Refrigerator 


Volume of refrigerator = 45 liter (42x42x26.5 cm^) 
Thickness of the glass wool insulation = 8 cm. 
Food item to be kept in refrigerator = 4 kg/h 
Food item always present in refrigerator = 15 kg 
Heat of respiration = 16 kJ/kg-day 
Specific heat of food items = 4.2 kJ/kg-deg.C 
Ice production = 12 kg/day 
Temperature of refrigerator = 2 deg. C 
Temperature of freezer = -10 deg. C 
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Appendix J 


Calculations of COP for 
Refrigerant R-12 

The p-h chart is shown in Figure (J.l)- From the property values of 
R- 12 , we have, 

enthalpy at point 1 = hi = 184.22 kJ/kg (J.l) 

entropy at point 1 = Si = 0.7060 kJ/kg-K (J-2) 

specific heat of R -12 vapourcp^ = 0.6 kJ/kg 
Considering a superheat of 10 °C from 1 to 1 ', we have, 
h[ = hi + lOcpg = 190.22 kJ/kg and (J.3) 

s; = Si -f Cpgi InTi/Ti = 0.7060 + 0.6 In 273.15/263.15 = 0.7284. (J.4) 

We assume that the efficiency of compressor is 0 . 8 . 

The ideal point after compression is 2 ' , while the actual point is 2 . 

Now, s'a = s'l = 0.7284 = S 3 + Cp^ InT^/T^ = 0.6863 + 0.6 lnT^/318.15 (J.5) 

which givesTg = 341.28 -K and 

/i' =hs + Cpgm - Ts) = 220.338 kJ/kg (J. 6 ) 

Ideal compressor work = h 2 “ = 30.118 kJ/kg (J.7) 


146 



Actual work= 30.118/.9 = 33.46 kJ/kg (J.8) 
Therefore, h 2 = 190.22 + 33.46 = 223.68 kJ /kg (J.9) 

Assuming 5°C subcooling, 

= — bcpi = 79.90—5 x 1.0= 74.94 kJ/kg (J.IO) 

Cooling effect= h[ -h^ = 115.32 kJ/kg (J.ll) 

COP = cooling effect/work done= 115.32/33.46 = 3.405 
Therefore we take COP = 3.5 


(J.12) 
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Appendix K 


Characteristic Function 


The process by which individuals from the universal set X are determined to be either 
members or non members of a set can be defined by a characteristic or discrimination 
function. For a given set A this function assigns a value {jLa{x) to every x€X such that 


fJLAix) 


1 if and only if a; G A 
0 if and only if x ^ A 


(K.l) 


Thus the function maps elements of the universal set to the set containing 0 and 1. 
This can be indicated by 


fJ,A • -X =>• {0, 1} 


iK.2) 


Fuzzy set 

The characteristic function of a crisp set assigns a value of either 1 or 0 to each in- 
dividual in the universal set. This function can be generalized such that the values 
assigned to the elements of the universal set fall within a specified range and indicate 
the membership grade of these elements in the set in question. Larger values denote 
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higher degrees of set membership. Such a function is called a membership function and 
the set defined by it is a fuzzy set. 

If Xi is an element of the universal set and yUi is its grade of membership in A, then 
A is written as 


A = /ii/rci +)U2/a:2+,.. .,+/Xn/2:n (A'.S) 

where the slash is employed to link the elements of the support with their grades of 
membership in A and the sign indicates that the listed pairs of elements and 
membership grades collectively form the definition of the set A. For the case in which 
a fuzzy set A is defined on a universal set that is finite and countable, we may write 

A = f:tXi/x, {KA) 

i=l 

Similarly, when X is an interval of real numbers, a fuzzy set A is often written in the 
form 


A = 



(K.5) 


Fuzzy union 

The union of fuzzy sets A and B is denoted by AUB which results in the maximum 
membership grades: 


A[jB = I [//A (a;) V 
where V is interpreted as “maximum” . 
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Fuzzy intersection 

The intersection of fuzzy sets A and B is denoted by AnB (A and B) which results in 
the minimum membership grades: 

^ n ^ A {K.7) 

where A is interpreted as “minimum” . 


Fuzzy product 

The product of A and B is denoted by AB which results in the ordinary product of 
two membership grades: 


AB = 


L 




{K.8) 


Convexity 

An a-cut of a fuzzy set A is a crisp set Aa that contains all the elements of the universal 
set X that have a membership grade in A greater than or equal to the specified value 
of a. This definition can be written as 

Aa = [xeX\fiA{x) > a] {K.9) 

A fuzzy set is convex if and only if each of its n-cuts is a convex set. Equivalently, 
we may say that a fuzzy set A is convex if and only if 


//^(Ar + (1 - A)s) > min /^yi(s)] for all r, se X and all Ae[0, 1] (if.lO) 
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Normality 

The fuzzy set A is called normal if the upper bound of its membership grades over X 
is unity, i.e. 


Sup xI^a{x) = 1 


{K.ll) 


Fuzzy numbers 

Fuzzy arithmetic is used for computations in which imprecise data exists in a simple or 
complex system. A convex and normalized fuzzy set defined on R (set of real numbers) 
whose membership function is piecewise continuous is called a fuzzy number. Thus, a 
fuzzy number can be thought of as an interval containing the real numbers to varying 
degrees. Figure K.l shows some examples of fuzzy numbers. 


Operations on fuzzy numbers 

(a) Fuzzy addition-addition of two fuzzy intervals at an a-cut in R is defined by 

Aa(+)B^ = K, aJ](+)(Sf , fif] = (af + 6f , of + ij] (Jf.W) 

(b) Fuzzy subtraction-subtraction of two fuzzy intervals at an a-cut in R is defined 
by 

= [<, <](-)[6r, bt] = K - m (K.n) 

(c) Fuzzy multiplication-multiplication of two fuzzy intervals at an o'-cut in R^'''^ is 

defined by 

= [a?, a^](x)[fe^ 6“] = [a? x ^2 x &?] (i^.l4) 
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X? 



b 


•D 



C9 

U 



Fig*K-1 Membership functions of (a) (o“y)o i 

(b) b j (c) n j(d) coefficient of friction 



(d) Fuzzy division-division of two fuzzy intervals at an a-cut in is defined by 


^a(-^)-Sa 




62 ) 0^2 




(ir.l5) 
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Appendix L 


Complex Search Method 


This method was developed by M. J. Box in 1965. The algorithm begins with a number 
of feasible points created at random. If a point is found to be infeasible, a new point 
is created using the previously-generated feasible points. Usually, the infeasible point 
is pushed towards the centroid of the previously-found feasible points. Once a set of 
feasible points is found, the worst point is reflected about the centroid of rest of the 
points tp find a new point. Depending on the feasibility and function value of the 
new point, the point is further modified or accepted. If the new point falls outside the 
variable boundaries, the point is modified to fall on the violated boundary. If the new 
point is infeasible, the point is retracted towards the feasible points. The worst point is 
replaced by this new feasible point and the algorithm continues for the next iteration. 

Algorithm 

Step 1 : Assume a bound in a reflection parameter a and termination 

parameters e, S. 

Step 2 : Generate an initial set of P (usually 2N) feasible points. For each point 
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(a) Sample N times to determine the point in the given bound. 

(b) If is infeasible, calculate x (centroid) of current set of points and reset 

until x^^ is feasible; 

Else if x^^^ is feasible, continue with (a) until P points are created. 

(c) Evaluate f{x^^^) forp = 0,1,2, , (P-1). 

Step 3 : Carry out the reflection step: 

(a) Select x^ such that 

f{x^) = maxf{x^'>) = Fmax- 

(b) Calculate the centroid x (of points except x^) and the new point 

x'^ =X + Oi{x - x^). 

(c) If is feasible and fix"^) > Fmax, retract half the distance to the centroid x. 

Continue until f{x^) < Fmax 

Else if is feasible and fix^) < Fmax, go to Step 5. 

Else if x'^ is infeasible, go to Step 4. 

Step 4 : Check for feasibility of the solution 

I 

(a) For all i, reset violated variable bounds: 

If xf^ < Xi^'^ set x'^ = xf'^ . 

If xf < xf^'^ set xf = xf^\ 
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(b) If the resulting x'^ is infeasible, retract half the distance to the centroid until x^ 


is feasible. Go to Step 3(c). 

Step 5 : Replace x^ by x^. Check for termination. 

(a) Calculate / = and x = jT,pX^\ 

(b) If -if <( and y'Epl|i“-i|P < S 

Terminate 

Else set fc = ^: + 1 and go to Step 3(a). 
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